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Using numerical simulations we investigate the properties of the dynamic phase transition that 
is encountered in the three-dimensional Ising model subjected to a periodically oscillating magnetic 
field. The values of the critical exponents are determined through finite-size scaling. Our results 
show that the studied non-equilibrium phase transition belongs to the universality class of the 
equilibrium three-dimensional Ising model. 
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I. INTRODUCTION 

In the last fifty years our understanding of equilibrium 
critical phenomena has developed to a point where well 
established results are available for a range of different 
situations. In particular, the origin of and the differ- 
ence between equilibrium universality classes is now well 
understood. Far less is known, however, on the proper- 
ties of non-equilibrium phase transitions taking place in 
interacting many-body systems that are far from equi- 
librium. As for the equilibrium situation, a continuous 
non-equilibrium phase transition is also characterized by 
a set of critical indices. Absorbing phase transitions 
encountered in reaction-diffusion systems^ and phase 
transitions found in driven diffusive systems^ are well 
studied examples of continuous phase transitions taking 
place far from equilibrium. Even though the values of the 
critical exponents have been determined in many cases, 
a classification of non-equilibrium phase transitions into 
non-equilibrium universality classes is far from complete. 

A different type of non-equilibrium criticality is en- 
countered when kinetic ferromagnets are subjected to a 
periodically oscillating magnetic field 5 ' 5 . When increas- 
ing the frequency of the field, a phase transition takes 
place between a dynamically disordered phase at low fre- 
quencies, where the ferromagnet is able to follow the 
changes of the field, and a dynamically ordered phase 
at high frequencies, where the magnetic system does not 
have time to adjust to the magnetic field before it changes 
its orientation. In the dynamically disordered phase the 
magnetization oscillates around zero, but in the dynam- 
ically ordered phase the magnetization keeps oscillating 
around a non- vanishing average value. This type of phase 
transition has been studied in a range of systems, with 
applications in chemical, physical, and materials systems 
as well as in ecology, and this both theoretically^— and 
experimentall y 14 i 15 . Most of these studies looked at gen- 
eral properties of the dynamical ordering, without trying 
to fully characterize the corresponding continuous tran- 
sition through the determination of the critical indices. 
It is only for the two-dimensional kinetic Ising model 
in a periodically oscillating field that the critical expo- 
nents have been measured^. Interestingly, even though 



the studied phase transition is a non-equilibrium transi- 
tion, the values of the exponents were found to coincide 
with the values of the universality class of the equilib- 
rium two-dimensional Ising model. This result agrees 
with an earlier conjecture^, based on the symmetry of 
the problem, as well as with an investigation of the time- 
dependent Ginzburg-Landau model with a periodically 
changing fielcUi From the symmetry argument given in 
Ref. [16] one expects also for other cases that this dy- 
namic phase transition falls into the universality class of 
the equilibrium counterpart, but this has not yet been 
confirmed. Interestingly, we showed in a recent studyi^ 
that the above mentioned symmetry argument breaks 
down when the dynamic phase transition takes place in 
systems with boundaries, yielding surface critical expo- 
nents at the non-equilibrium phase transition that do not 
coincide with the corresponding equilibrium exponents. 

In this short communication we present evidence that 
also for other bulk systems the dynamic phase transition 
falls into the corresponding equilibrium universality class. 
This is done through a study of the critical properties of 
the kinetic three-dimensional Ising model in a square- 
wave field. 

In the next Section we briefly discuss the three- 
dimensional kinetic Ising model and the quantities used 
to characterize the critical properties at the dynamic 
phase transition. Section III contains the results of our 
numerical study. We conclude in Section IV. 



II. MODEL 

The Hamiltonian of the three-dimensional kinetic Ising 
model is given by the expression 



(1) 
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where S x = ±1 is the usual Ising spin located at site x. 
The first sum in that expression runs over bonds con- 
necting nearest neighbor spins. J > is a ferromagnetic 
coupling and H(t) is a spatially uniform periodically os- 
cillating magnetic field. Following Q we use a square- 
wave field of amplitude Hq and half-period ti/2- For the 
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data presented in this paper we choose temperature and 
magnetic field strength so that the system is in the mul- 
tidroplet regime: T = 0.8 T c and H = 0.4 J, with the 
critical temperature T c — 4.5115J//cb (ks is the Boltz- 
mann constant). 
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FIG. 1: (Color online) Time-dependent magnetization of the 
three-dimensional kinetic Ising model in a square-wave field 
of half-period (a) t 1/2 = 100, (b) t 1/2 = 60, and (c) t 1/2 = 20. 
Starting with the fully magnetized sample, the first five peri- 
ods are shown for each case. The red dashed line is H(t)/Ho 
where Ho is the amplitude of the magnetic field. Whereas for 
case (a) we are in the dynamically ordered phase, for case (c) 
the system is dynamically disordered. Case (b) is close to the 
dynamical phase transition. The parameters are T = 0.8 T c 
and Ho = 0.4 J, the system contains 96 s spins. The unit of 
time is one Monte Carlo Step. 

Let us pause briefly in order to summarize the mecha- 
nism underlying the dynamical ordering observed in the 
kinetic ferromagnets, see Fig. [TJ Assuming the magne- 
tization being aligned with the external field, a revers- 
ing of the direction of the field results in a metastablc 
state from which the system tries to escape by nucleat- 
ing droplets that are pointing in the same direction as the 
field. Whereas in the dynamically disordered phase, the 
ferromagnet is able to reverse its magnetization before 
the field is changed again, see Fig. Q]}, in the dynamically 
ordered phase the system is still in the metastable state 
when the field direction switches back, thereby yielding 
a time dependent magnetization that oscillates around 
a finite value, as shown in Fig. [T];. This competition 
between the magnetic field and the metastable state is 
captured by the variable 



where (r) is the metastable lifetime. For the kinetic Ising 
model plays the same role as that played by temper- 
ature in the equilibrium system. Changing t±/2 changes 
the value of 0. Consequently, there is a critical value 
C where the dynamic phase transition between the dy- 
namically ordered and the dynamically disordered phases 
takes place, see Fig. [T}d. 

In order to determine the metastablc lifetime we con- 
sider an equilibrated system in a field of constant strength 
Ho. Reversing the sign of the field renders the system 
metastable. The metastable lifetime is then obtained as 
the average first-passage time to zero magnetization, see 
Fig. [21 For our system parameters, T = 0.8 T c and 
Hq = 0.4 J, we obtain that (r) = 47.05, where the time 
is measured in Monte Carlo Step, with one Step corre- 
sponding to updating on average every spin once. 
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FIG. 2: (Color online) Determination of the metastable life- 
time. After preparing the system in a fully ordered state a 
constant magnetic field is applied that reverses the magneti- 
zation. The metastable lifetime is defined as the first-passage 
time to zero magnetization. The parameters are T = 0.8 T c 
and Ho = 0.4 J, the data shown have been obtained for 
a system composed of 96 3 sites. The metastable lifetime 
(r) = 47.05 is independent of the system size. The unit of 
time is one Monte Carlo Step. 

In order to elucidate the critical properties of the dy- 
namic phase transition we study a range of different 
quantities. The main quantity of interest is the period- 
averaged magnetization 

Q = J— I m(t)dt (3) 

**l/2 J 

where the integration is performed over one period of 
the oscillating field. The time dependent magnetization 
is given by 

X 

We thereby consider cubic systems with periodic bound- 
ary conditions and linear extend L, so that N = L 3 . The 
order parameter is then obtained after averaging 
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both over time (i.e. an average over many periods) and 
over different realizations of the noise. 

Some quantities of interest are directly obtained from 
the order parameter. Thus we calculate the Binder cu- 
mulant 



17 = 1- 



(Q 4 ) 
HQ 2 ) 2 



(5) 



that we use to determine the critical point, see Section 
III. The susceptibility describing the fluctuations of the 
order parameter is given by 



X Q =iV«Q 2 )-<|Q|> 2 ) . 
We also calculated the period-averaged energy 

(x>y) 

and its fluctuations 

X E = N((E 2 ) - (E) 2 ) . 
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Close to the critical point, the order parameter and 
the response functions x Q an d X E display an algebraic 
behavior: 
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with the critical exponents /3, 7, and a. 



III. RESULTS 

As a necessary prerequisite for a reliable determination 
of the critical exponents we need to know the location of 
the critical point with high precision. As usual for con- 
tinuous phase transitions the various quantities show a 
characteristic behavior in finite systems when approach- 
ing the critical point. This is illustrated in Fig. [3] for the 
Binder cumulant, the order parameter, as well as for the 
fluctuations of the order parameter. Even though one can 
get an estimate for the critical value C by analyzing the 
shift of the maximum of the susceptibility when changing 
the system size, the most reliable way to determine the 
critical point is to study the Binder cumulant. As seen in 
Fig. [3Ji, for the kinetic three-dimensional Ising model the 
data sets obtained for different system sizes all cross at 
a common value = 1.285. This value being compatible 
with the finite-size shifts of the position of the maxima 
in and x E ■> we take this value as the critical value C 
where the dynamic phase transition takes place. 

In order to determine the values of the critical expo- 
nents, we perform a finite-size scaling analysis^, similar 
to what has been done previously for the two-dimensional 
kinetic Ising model^. Close to the critical point systems 
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FIG. 3: (Color online) (a) Binder cumulant, (b) order param- 
eter, and (c) susceptibility as a function of for systems with 
different linear extend L. The cumulants for the different sys- 
tem sizes cross at Q c = 1.285. In addition, the susceptibility 
displays in the finite systems a maximum close to that value 
of G. The data result from averaging over typically 100000 
periods. 



of different sizes should display the following scaling be- 
havior: 



{\Q\) = L-PI V T ± {9L 1 '") 
Li' v g±{6L 1 ' v ) 



,Q - 



X 



= L a l v X ± {dL 1 l v ) 



(12) 

(13) 
(14) 



where 6 = 



e-e c 



is the reduced control parameter and v 



is the critical exponent describing the divergence of the 
correlation length when approaching the critical point. 
F±, Q±, and T± are scaling functions, where the sign ± 
correspond to ^ C . Comparison of data obtained 
at the critical point = C should yield for all three 
quantities an algebraic dependence on the system size, 
which allows to directly determine the critical exponents 
from the slopes in a double logarithmic plot. 
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FIG. 4: (a) Order parameter (\Q\), (b) order parameter sus- 
ceptibility 'j an d (c) energy fluctuations \ E a t the critical 
point 9 = 6c as a function of system size. The slopes in 
a double logarithmic plot give the values of the critical ex- 
ponents P/v, 'y/v, and a/ v. The data have been obtained 
by averaging over at least 200000 periods. System sizes are 
measured in multiples of the lattice constant. 

As shown in Fig. @] all three quantities indeed depend 
algebraically on the system size. From the slopes we 



4 



obtain the values lijv = 0.51(2), j/u = 1.96(4), and 
a/v = 0.18(3) (for the specific heat we did not take into 
account the largest system size, due to the rather poor 
quality of the data). Comparing these values with the 
well known equilibrium values for the three-dimensional 
Ising model [3 = 0.33, 7 = 1.24, a = 0.108, and v = 0.63, 
we see that the set of exponents we obtain for the dy- 
namic phase transition of the kinetic Ising model in- 
deed agree with the equilibrium values j3/v = 0.52, 
7/V = 1.97, and a/v — 0.17. This confirms that for 
the problem at hand the symmetry argument given in 
Ref. [lj| is also valid in three space dimensions. 

IV. CONCLUSION 

Classifying and understanding universality classes of 
non-equilibrium phase transitions remains an important 
task. For some types of systems, as for example absorb- 
ing phase transitions, field theoretical methods can be 
used successfully, which then allow not only to calculate 
critical quantities but also to understand what features 
of a given system are universal. Still, the classification of 
non-equilibrium universality classes is far from complete, 
and much work remains to be done before a more com- 
plete understanding of the origin of universality far from 



equilibrium is achieved. 

The dynamic phase transition encountered in magnetic 
systems in an oscillating field constitutes an interesting 
case. Indeed, it was shown that for the two-dimensional 
kinetic Ising model the critical exponents at this non- 
equilibrium phase transition coincide with those of the 
equilibrium Ising models. This result is supported by a 
symmetry argument^. On the other hand, the presence 
of surfaces results in non-equilibrium surface universality 
classes that differ from those encountered in the equilib- 
rium semi-infinite Ising models. 

In this work we have shown that also for the dy- 
namic phase transition in the three-dimensional kinetic 
Ising model the same critical exponents as for the three- 
dimensional equilibrium model are obtained. This lends 
further support to the symmetry argument given in [l6| . 
It also stresses the importance of symmetries in establish- 
ing universalities in non-equilibrium critical phenomena. 
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